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Abstract 

Two choices are possible for the refractive index of a linear, homogeneous, isotropic, active, 
dielectric material. Either of the choices is adequate for obtaining frequency-domain solutions for 
(i) scattering by slabs, spheres, and other objects of bounded extent; (ii) guided-wave propaga- 
tion in homogeneously filled, cross-sectionally uniform, straight waveguide sections with perfectly 
conducting walls; and (iii) image formation due to flat lenses. The correct choice does matter 
for the half-space problem, but that problem is not realistic. 



1 Introduction 

A century and a half after the unification of light with electricity and magnetism, certain aspects of 
electromagnetic wave propagation in a linear, homogeneous, isotropic, dielectric material continuum 
remain unsettled. One of those aspects is as follows. If e r (uj) = e' r (ui) + £e"(u>) denotes the relative 
permittivity — where e' r 6 R, e" 6 K, and an exp(— iut) time-dependence is implicit — then the 
refractive index n(u>) satisfies the dispersion relation n 2 (w) = e r (uj). Is the refractive index equal to 
71+ or 7i- = — n+ , where Im(n+) > 0? 

Suppose that a material is passive at a specific frequency Q, i.e., e"(Q) > 0. The general consensus 
is that the refractive index is equal to 7i+(f2). If the material is active at the frequency (i.e., 
e"(S7) < 0), however, disagreements on the correct choice of the refractive index have surfaced. We 
are exclusively concerned with linear, homogeneous, isotropic, active, dielectric materials in this 
communication, subject to the additional stipulations that the constitutive properties are spatially 
local and unalterable by the passage of an electromagnetic signal. 

Many reasons have been provided for the refractive index of such a material to be 7i_(f2) [l]-[7] 
and n+ (17) [8]— [11]. A direct time-domain solution of the Maxwell equations has recently upheld 
the possibility of n+(fi) as the refractive index [12]. Even more importantly, the adoption of the 
time-domain procedure implies that the refractive index cannot be determined at a single frequency 
by frequency-domain analyses, if the material is active at that frequency. By virtue of this argument, 
it is possible for two different active materials to have the same permittivity but opposite refractive 
indexes at a certain frequency [13]. 

Thus, in order to choose the refractive index of an active material at a specific frequency, the 
solution of a time-domain problem is needed with an incident signal of sufficient bandwidth to cover 



the frequency of interest (though it appears that a frequency-domain algorithm may sometimes 
suffice to determine the variation of n(u) vs. lo over a sufficiently wide w-range [13]-[15]). The time- 
domain problem should be such that the correct choice of the refractive index in the corresponding 
frequency-domain problem must be consequential — e.g., the reflection of a signal by a half-space 
filled with a material continuum [12]. This naturally leads to the question: when is knowledge of the 
correct choice of the refractive index essential for the solution of a frequency-domain problem? 

In the following section, several commonly tackled frequency-domain boundary-value problems 
are discussed, in order to answer the foregoing question. The selected problems include scattering 
by objects of bounded extent; guided-wave propagation in homogeneously filled straight waveguide 
sections with perfectly conducting walls and uniform cross-section; and lensing by slabs. 



2 Frequency— domain Boundary— value Problems 



2.1 Planewave response of a slab 

Let us begin by considering the electromagnetic response of the region < z < L filled with an active 
material. The regions z < and z > L are vacuous. 

Because the dyadic Green function for free space can be expanded in terms of an angular spectrum 
of plane waves, we can consider the response of the slab to a plane wave without significant loss of 
generality. An incident plane wave can be represented by 



E inc (r,w) 
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exp [i (kx + aoz)] , z < , 
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where a s and a p are the amplitudes of the s- and the p-polarized components, respectively; fco is the 
free-space wavenumber; and a = +\Ao — k2 wrt h K S The reflected and the transmitted plane 
waves are represented by 



and 
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respectively, while the electric field phasor inside the slab region is 



exp {i [kx + q:q(z — L)]} , z > L . 
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exp [i (kx — az)] , < z < L . 
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where a = +y/k^n 2 — k 2 . Standard techniques yield the amplitudes 

Tp = ~a p (e 2taL - l)(a 2 - n 4 a 2 )/A p , t p = -4n 2 a p e 2laL aao/ A p 
c p = —2a p na (a + n 2 ao)/A p , d p = — 2a p nao(a — n 2 ao)/A p 

r s = -a s (e 2laL - l)(a 2 - a 2 )/A s , t s - -4a s e 2laL aa / A s 
c s = -2a s a (a + a )/A s , d s = 2a s a (a - a )/A s 



(5) 
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where 



A p = (e 2iaL - l)(a 2 + n 4 ag) - 2n 2 {e 2iaL + l)aa 
A s = (e 2laL - l)(a 2 + a 2 ) - 2{e 2laL + l)aa 




The reflected, transmitted, and the internal electric field phasors are not affected by the change 
n — > — n; neither are the corresponding magnetic field phasors. Therefore, the fields everywhere are 
unaffected whether n = n + or n = n_ is chosen. The same conclusion is obtained if the plane z = L 
is assumed to be perfectly conducting. 



2.2 Planewave response of a sphere 



Let us next consider the spherical region r < a occupied by an active material, whereas the region 
r > a is vacuous. Without loss of generality, we take the sphere to be illuminated by a linearly 
polarized plane wave traveling along the +z axis. As is commonplace, the incident plane wave is 
represented in terms of vector spherical harmonics, M.amui'w) and N^rr^w) [16], as 



E mc (r,c) = AJ2 ^StH Kit (M - iN&for) 



(7) 



where A is the amplitude. The scattered and the internal electric field phasors are also expressible 
in terms of vector spherical harmonics; thus [16], 
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ji/(0 an d h\}\^) are the spherical Bessel function and the spherical Hankel function of the first kind, 



respectively; and 



(14) 



Expressions for the corresponding magnetic field phasors can be derived using the Faraday equation. 
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Let the change n — > — n be effected in Eqs. (7)-(13). Because = (— and^i* (— £) = 

(— l)^^! 1 ' (^), it follows that both a v and i>„ are not affected, and neither are the scattered electric 
and magnetic field phasors. Although c v — ► (— 1)" c v and — ► (— l) v+1 d Vl the internal electric and 
magnetic field phasors are not affected because Mct2j„(— w) = (-l)"Mm„(w) and N^i>(— w ) = 
(-l)' y+1 N^(w). Therefore, again, the fields everywhere do not depend on the choice of the 
refractive index of the active material. We have numerically verified that the same conclusion holds 
for a dielectric sphere with a perfectly conducting, concentric, spherical core. 



2.3 More general scattering problems 

The conclusion garnered in Sec. 2.2 is valid in a far more general situation. Let all space be divided 
into three mutually disjoint regions Vq, Vj, and V e as follows. Whereas both Vj and V e are bounded 
in extent, the vacuous region Vo extends to infinity in all directions. The region V e is filled with an 
active material, whereas the sources of the electromagnetic field reside wholly in Vj. 

The solution of the frequency-domain Maxwell curl equations everywhere can be stated as [17] 
E(r, uj) = iuj^i a / Q{k r, k r') • J so (r', w) d 3 r' 

JVj 

+ k 2 [ [n 2 (w)-l]^(fc r,fc r') 'E(r»rfV, reVjU^UVb, (15) 
Jv e 

where jj, is the free-space permeability, G(fcor, kor') is the dyadic Green function for free space, and 
J so (r, uo) is the source electric current density phasor. An expression for the magnetic field phasor can 
be derived from Eq. (15) and the Faraday equation. Once again, the electric and the magnetic field 
phasors everywhere are not affected by the change n — ► —n. The same conclusion is obtained if the 
active material is nonhomogeneous, for which case n 2 (ui) must be replaced by n 2 (r',uj) in Eq. (15). 



2.4 Homogeneously filled waveguides of uniform cross-section 



A commonplace problem in electromagnetism is the propagation of waves in straight waveguides of 
uniform cross-section and perfectly conducting walls. Let us begin with single-conductor waveguides. 
An example is the rectangular waveguide whose walls are formed by the planes x = 0, x = a, y = 0, 
and y = b. Let the propagation direction be along the +z axis. The electric field phasor in this 
waveguide may be represented as the modal sum [18] 



E 



(r,w) = A m E) cx p(>9 z ) 7^ ~y cos sin * 
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wherein the combination p = q = is not allowed; A P ^ E ^ and B^ M ^ are modal coefficients for 
transverse-electric and transverse-magnetic modes, respectively; and 7^ = A pq — k$ n 2 with A pq — 
(pn/a) 2 + (qir/b) 2 . A similar expression for magnetic field phasor can be derived using the Faraday 
equation. The value of the propagation constant 7 pg of the (pq) th mode is chosen to ensure that the 
integral of the modal time-averaged Poynting vector over any xy plane is coparallel with z. Neither 
the electric nor the magnetic field phasor is affected by the change n — > —n. The same conclusion 
holds true for parallel-plate waveguides (b — > 00) and circular waveguides [18]. 

More generally, in a cylindrical coordinate system (u, v, z), let E(r, u) = [e t (u, v, u) + e z (u, v, u>)z] 
exp(7z) and H(r, u>) = [h t (u, v, u>) + h z (u, v, u)z] exp(7z), where z • e t = and z • h t = 0. Then, the 
source-free Maxwell curl equations yield the transverse components of the field phasors in terms of 
the longitudinal components as [18] 



e*(«,v,w) - 
Mu,v,uj) = 



(grad t e 2 -»zx (grad t /i z ) 



(grad^ + z&B- z x (grad t e 2 ) 



(17) 



where grad t = grad — z {d/dz}. The longitudinal field components satisfy the Hclmholtz equation 
as follows: 

For guided-wave propagation, 7 2 emerges from a dispersion equation derived from enforcing the 
usual boundary conditions on the perfectly conducting walls. Clearly, to obtain modal solutions for 
e z and h z , the refractive index is not needed by itself, but its square is needed. Therefore, the fields 
arc unaffected whether n = n + or n = n_ is chosen. 

In a two-conductor waveguide, in addition to transverse-electric and transverse-magnetic modes, 
a transverse-electromagnetic mode can also propagate along the +z axis. This latter type of mode 
is like a plane wave, and its description does require knowledge of the correct refractive index. 
The number of transverse-electromagnetic modes increases with the number of conductors [18]. As 
waveguide sections of finite length are used in practice, modal propagation along the — z and +z 
directions must be considered simultaneously [19]. The modal coefficients of the counterpropagating 
modes in a waveguide section would adjust according to the boundary conditions enforced at both 
ends of the section, just as in Sec. 2.1 for a slab; thus, the correct choice of the refractive index need 
not be critical. 



2.5 Flat lenses 

A major motivation for current research on isotropic dielectric-magnetic materials that display neg- 
ative phase velocity (i.e., Re(n) < 0) is their potential to form flat lenses. Ideally, a flat lens is a slab 
made of nondissipative and nondispersive material with e r = — 1 and relative permeability /j, r = — 1 
(i.e., anti-vacuum). When this lens is sandwiched between vacuous half-spaces, an image of a source 
placed on one side of the lens is formed on the other side of lens and another image is formed within 
the lens itself. Both propagating and evanescent modes are brought to a focus at the image points 
in this ideal scenario [20]. A rigorous frequency-domain analysis reveals that the creation of images 
is independent of the choice of the refractive index of the lens material; furthermore, the focusing 
properties are seriously compromised if e r = /i r 7^ —1 [21]. The imaging capabilities are further com- 
promised by dissipation and/or dispersion, and the extent to which the slab demonstrates paraxial 
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focusing or channeling of the illuminating field is independent of whether the slab's refractive index is 
chosen to lie in the upper or lower half-space in the complex plane [21]. From these results, it follows 
that the nonmagnetic, active dielectric material of our interest here is unsuitable for flat lenses. 



3 Discussion 

There remains the problem of reflection of a plane wave from a half-space filled with an active 
material [13, 15]. As has been shown by reference to the time-domain solution [12], only one of 
the two choices of the refractive index yields the correct frequency-domain reflectance; that is also 
the conclusion from an analysis with the Laplace transform [15]. However, a time domain solution 
can never access an entire half-space; thus, the removal of the ambiguity through a comparative 
study may not be considered by some researchers to be totally compelling. Such a comparative 
study may not be possible for many active materials anyway, because the fields will continue to grow, 
invalidate the assumption of linearity, and eventually counter the assumption of temporal stationarity 
of the constitutive properties. Furthermore, whereas when solving a frequency-domain problem, a 
half-space can be considered to be an adequate simplification of a sufficiently thick slab when the 
material is dissipative, the process of taking the limit L — > oo in Sec. 2.1 — for an active material 
[10, 11] — has unphysical consequences in general [15]. This is because the back face of an active 
slab can receive a sufficiently strong field which it can reflect back to reach the front face, but a 
half-space has no back face. 

To conclude, we have shown that, for the most common frequency-domain problems handled 
by electromagnetics researchers, knowledge of the correct choice of the refractive index of a linear, 
homogeneous, isotropic, active, dielectric material is inessential. Either of the two choices would 
serve adequately. The major exception we encountered is the half-space problem, but that problem 
is significantly unrealistic. 
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